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We consider a host-parasite model for a population of cells that can be 
of two types, A or B, and exhibits unilateral reproduction: while a B-cell 
always splits into two cells of the same type, the two daughter cells of an 
A-cell can be of any type. The random mechanism that describes how 
parasites within a cell multiply and are then shared into the daughter cells 
is allowed to depend on the hosting mother cell as well as its daughter cells. 
Focusing on the subpopulation of A-cells and its parasites, our model differs 
from the single-type model recently studied by Bansaye [5] in that the 
sharing mechanism may be biased towards one of the two types. Our main 
results are concerned with the nonextinctive case and provide information 
on the behavior, asn-> oo, of the number A-parasites in generation n and 
the relative proportion of A- and B-cells in this generation which host a 
given number of parasites. As in [5], proofs will make use of a so-called 
random cell line which, when conditioned to be of type A, behaves like a 
branching process in random environment. 

1 Introduction 

The reciprocal adaptive genetic change of two antagonists (e.g. different species or 
genes) through reciprocal selective pressures is known as host-parasite coevolution. 
It may be observed even in real-time under both, field and laboratory conditions, if 
reciprocal adaptations take place rapidly and generation times are short. For more 
information see e.g. QTJ[T2]. 

The present work studies a host-parasite branching model with two types of cells (the 
hosts), here called A and B, and proliferating parasites colonizing the cells. Adopting 
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a genealogical perspective, we are interested in the evolution of certain characteristics 
over generations and under the following assumptions on the reproductive behavior of 
cells and parasites. All cells behave independently and split into two daughter cells 
after one unit of time. The types of the daughter cells of a type-A cell are chosen 
in accordance with a random mechanism which is the same for all mother cells of 
this type whereas both daughter cells of a type-B cell are again of type B. Parasites 
within a cell multiply in an iid manner to produce a random number of offspring the 
distribution of which may depend on the type of this cell as well as on those of its 
daughter cells. The same holds true for the random mechanism by which the offspring 
is shared into these daughter cells. 

The described model grew out of a discussion with biologists in an attempt to provide 
a first very simple setup that allows to study coevolutionary adaptations, here due to 
the presence of two different cell types. It may also be viewed as a simple multi-type 
extension of a model studied by Bansaye [5] which in turn forms a discrete-time 
version of a model introduced by Kimmel [TU] . Bansaye himself extended his results 
in [6] by allowing immigration and random environments, the latter meaning that each 
cell chooses the reproduction law for the parasites it hosts in an iid manner. Let us 
further mention related recent work by GUYON [5] who studied another discrete-time 
model with asymmetric sharing and obtained limit theorems under ergodic hypotheses 
which, however, exclude an extinction-explosion principle for the parasites which is 
valid in our model. 

We continue with the introduction of some necessary notation which is similar to 
the one in [5]. Making the usual assumption of starting from one ancestor cell, denoted 
as 0, we put Go := {0}, G n := {0, l} n for n > 1, and let 

T := (J G„ with G„ := {0, 1}™ 

be the binary Ulam-Harris tree rooted at which provides the label set of all cells in 
the considered population. Plainly, G n contains the labels of all cells of generation n. 
For any cell »gT, let T v £ {A, B} denote its type and Z v the number of parasites it 
contains. Unless stated otherwise, the ancestor cell is assumed to be of type A and to 
contain one parasite, i.e. 

T = A and Z = 1. (SAl) 
Then, for t e {A, B} and n > 0, define 

G„(t) := {v e G„ : T v = t} and G^(t) := {v e G„(t) : Z v > 0} 

as the sets of type-t cells and type-t contaminated cells in generation n, respectively. 
The set of all contaminated cells in generation n are denoted G* , thus G* = G* (A) U 

g;(b). 

As common, we write V\...v n for v = (v\, ...,v n ) £ G n , uv for the concatenation of 
u, v G T, i.e. 

uv = u\..., u m v\...v n if u = u\...u m and v — v\...v ni 
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and v\k for the ancestor of v = vi...v n in generation k < n, thus v\k = v\,...,Vk- 
Finally, if v\k = u for some k and m/d, we write u < v. 

The process (T v ) v& j is a Markov process indexed by the tree T as defined in [7]. It 
has transition probabilities 

P(T vQ = x,T vl =y\T v = A) = p xy , (x,y) G {(A, A), (A, B), (B, B)}, 

P(T„ - B,T vl = B\T V = B) = 1, 

and we denote by 

Po ■= Paa + Pab = 1 - Pbb and pi := p A A 

the probabilities that the first and the second daughter cell are of type A, respectively. 
In order to rule out total segregation of type- A and type-B cells, which would just lead 
back to the model studied in [5_ , it will be assumed throughout that 

Paa < I- (SA2) 

The sequence (#G„(A))„> obviously forming a Galton- Watson branching process 
with one ancestor (as T = A) and mean 

v := po + Pi = 2paa + Pab = 1 + (paa - Pbb) < 2, 

it is a standard fact that (see e.g. [4]) 

#G„(A) -> a.s. iff paa < Pbb and p A B < 1. 

To describe the multiplication of parasites, let Z v denote the number of parasites in 
cell v and, for t g {A, B}, s g {AA, AB, BB}, let 

(4°](t,s),X«(t,s)) , t G {A, B}, s G {AA, AB, BB} 

be independent families of iid Ng-valued random vectors with respective generic copies 
(X(°)(t, s), X^^(t, s)). If v is of type t and their daughter cells are of type x and y, then 
xj?' v (t,xy) gives the offspring number of the k th parasite in cell v that is shared into 
the daughter cell v i of v. Since type-B cells can only produce daughter cells of the same 
type, we will write (xg(B), xg(B)) as shorthand for (^ ](B, BB), xg(B, BB)). To 
avoid trivialities, it is always assumed hereafter that 

P(^ (0) (A,AA) < 1, X (1) (A,AA) < l) < 1 (SA3) 

and 

P (x^(B) < 1, lW(B) < l) < 1. (SA4) 

Next, observe that 

te{A,B} se{AA,AB,BB} k=l 
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We put /ij,t(s) := EJ«(t,s) for i S {0,1} and t,s as before, write /U^b as shorthand 
for /^^(BB) and assume throughout that /Zj,t(s) are finite and 

Mo.aCAA), a*i,a(AA), E (# G*(B)) > 0, mo,b, /*i,b > 0. (SA5) 

The total number of parasites in cells of type t G {A, B} at generation n is denoted 

by 

Z n {t) := y ' 

»6G»(t) 

and we put Z„ := Z„(A)+Z„(B), plainly the total number of all parasites at generation 
n. Both, (Z n ) n >o and (2„(A))„>o, are transient Markov chains with absorbing state 
and satisfy the extinction-explosion principle (see Section 1.5 in [4. for a standard 
argument), i.e. 

P(Z„ -> 0) + P(Z„ ->• oo) = 1 and P(Z„(A) -> 0) + P(Z„(A) ->• oo) = 1. 
The extinction events are defined as 

Ext := {Z„ 0} and Ext(t) := {Z„(t) -»()}, t£ {A, B}, 

their complements by Surv and Surv(t), respectively. 

As in [5] , we are interested in the statistical properties of an infinite random cell line, 
picked however from those lines consisting of A-cells only. This leads to a so-called 
random A-cell line. Since B-cells produce only daughter cells of the same type, the 
properties of a random B-cell line may be deduced from the afore-mentioned work and 
are therefore not studied hereafter. 

For the definition of a random A-cell line, a little more care than in [5] is needed 
because cells occur in two types and parasitic reproduction may depend on the types of 
the host and both its daughter cells. On the other hand, we will show in Section [5] that 
a random A-cell line still behaves like a branching process in iid random environment 
(BPRE) which has been a fundamental observation in 5 for a random cell line in the 
single-type situation. 

Let U — (U n ) 7U zfi be an iid sequence of symmetric Bernoulli variables independent 
of the parasitic evolution and put V n :— U\...U n . Then 

=: V -)• Vi ->• V 2 ... -)• V n ... 

provides us with a random cell line in the binary Ulam-Haris tree, and we denote by 

T[ n ] = T Vn and Z[ n ] = Z Vn n>0, 

the cell types and the number of parasites along that random cell line. A random 
A-cell line up to generation n is obtained when T[ n ] = A, for then T[ fe i = A for any 
k — 0, ...,n — las well. As will be shown in Prop.[U the conditional law of (^[o], Z\n\) 
given Tr n i = A, i.e., given an A-cell line up to generation n is picked at random, equals 
the law of a certain BPRE (2'fe(A))fc>o up to generation n, for each neN. It should 
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be clear that this cannot be generally true for the unconditional law of {Zm\, -^[n]), 
due to the multi-type structure of the cell population. 

Aiming at a study of host-parasite coevolution in the framework of a multitype host 
population, our model may be viewed as the simplest possible alternative. There are 
only two types of host cells and reproduction is unilateral in the sense that cells of 
type A may give birth to both, A- and B-cells, but those of type B will never produce 
cells of the opposite type. The basic idea behind this restriction is that of irreversible 
mutations that generate new types of cells but never lead back to already existing 
ones. Observe that our setup could readily be generalized without changing much the 
mathematical structure by allowing the occurrence of further irreversible mutations 
from cells of type B to cells of type C, and so on. 

The rest of this paper is organized as follows. We focus on the case of non-extinction 
of contaminated A-cells, that is P(Ext(A)) < 1. Basic results on Z n (A), Zu, #G*(A) 
and # G* including the afore-mentioned one will be shown in Section [2] and be partly 
instrumental for the proofs of our results on the asymptotic behavior of the relative 
proportion of contaminated cells with k parasites within the population of all contami- 
nated cells. These results are stated in Section [3] and proved in Section 0J A glossary of 
the most important notation used throughout may be found at the end of this article. 

2 Basic Results 

We begin with a number of basic properties of and results about the quantities G* (A) , 
G*, Z n (A) and Z [n] . 

2.1 The random A-cell line and its associated sequence (Z[ n ])n>o 

In [5] , a random cell line was obtained by simply picking a random path in the infinite 
binary Ulam-Harris tree representing the cell population. Due to the multi-type struc- 
ture here, we must proceed in a different manner when restricting to a specific cell 
type, here type A. In order to study the properties of a "typical" A-cell in generation 
n for large n, i.e., an A-cell picked at random from this generation, a convenient (but 
not the only) way is to first pick at random a cell line up to generation n from the 
full height n binary tree as in [5] and then to condition upon the event that the cell 
picked at generation n is of type A. This naturally leads to a random A-cell line up 
to generation n, for A-cells can only stem from cells of the same type. Then looking 
at the conditional distribution of the associated parasitic random vector (Zmu Z[ n ]) 
leads to a BPRE not depending on n and thus to an analogous situation as in [5] ■ The 
precise result is stated next. 

Proposition 1. Let (Z n (A)) n >Q be a BPRE with one ancestor and iid environmental 
sequence (A„)„>i taking values in {£(X<W(A, AA)), C(X^ (A, AA)), C(X^ (A, AB))} 
such that 

p(Aj = £pT (0) (A,AB))) = ^ and P (a x = C(X^(A, AA))) = 
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for i G {0, 1}. Then the conditional law of (Z^, given Tj ra ] = A equals the law 
of (Z (A), Z n (A)), for each n > 0. 

Proof. We use induction over n and begin by noting that nothing has to be shown if 
n = 0. For n > 1 and (z , 2„) G Nq +1 , we introduce the notation 

C 20 ,..., Zn := {{Z [0] ,...,Z [n] ) = (z ,...,z n )} and C* 0i ... iIb := C ZOt ... tZn n {T [n] = A} 

and note that 

P (T W = A) = 2""E ( ]T 1{t„=a } ) = g)" , 

\vEG n / 

for each n G N, in particluar 

P(T W = A IT^!] = A) = p(T[n _ i]=A) = 2- 

Assuming the assertion holds for n — 1 (inductive hypothesis), thus 

P(C Z0 ,..., Zn _ 1 |T [n _i ] =A)=P(Z (A) = «o,-,^n-i(A)=^ n -i) 
for any (z n , 2 n _i) S NJf, we infer with the help of the Markov property that 

P ((■%>], — ! Z[ n ]) = (z , Z n )\T[ n ] = A) 
P(<ff ,...,«„) 

~ P(T [n] =A) 

P(T [n -i] = A) 



= P(a„,..,z„_ 1 |T [n _ 1] =A) P(Z [n] =0„,T W =A|C, A 0i ... iZn _ 1 ) 
= F(Z (A)=z ,...,Z n _ 1 {A)=z n - 1 



P(T[n] = A) 
P(Z [1] = Zn,^!] = A|Z [0 ] = Z»„i,T [0 ] = A) 

P(T [n] = A \T [n _q = A) 
= F(Z (A)=z ,...,Z n _ 1 {A)=z n - 1 ) 

%feK )(A ' AB) )" _1 (^})+ E ^(p xW ^)) m "- 1 

= F(Z (A)=z ,...,Z n _ 1 (A)=z n - 1 )F(Z [n] (A) = z n \Z [n _ 1] (A) = z n - 1 ) 

= F(Z (A)=z ,...,Z n (A) = z n ). 

This proves the assertion. □ 

The connection between the distribution of Z n (A) and the expected number of A- 
cells in generation n with a specified number of parasites is stated in the next result. 
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Proposition 2. For all n eN and k 6 No, 

P (Z n (A) = fc) = i/-" E (#{« e G n (A) : Z„ = k}) , (1) 

m particular 

¥(Z n (A) >0) = 2/- n E#G;(A). (2) 
Proof. For all n, fc G N, we find that 

E(#{«eG„(A) :Z„ = A}) = ^ P(Z„ = k, T v = A) 

»eG„ 

= 2 n F(Z [n] =k,T [n] =A) 

- 2"P(T W = A)P (Z [n] = k\T [n] = A) 

= i/ n P (Z [n] = k\T [n] = A) 

= u n ¥(Z n (A) = k), 

and this proves the result. □ 

For neN and s E [0, 1], let 

/„(s|A) := E(,s z "( A '|A) and f n (s) := Es z »^ = E/ n (s|A) 

denote the quenched and annealed generating function of Z n (A), respectively, where 
A := (A n )„>i. Then the theory of BPRE (see P1I51I51[T2] for more details) provides us 
with the following facts: For each neN, 

/ n (s|A) = g Al o ... og An (s), g x (s) := E(«*W|Ai = A) = £ X n s n 

n>0 

for any distribution A = (A n ) n >o on No- Moreover, the g An are iid with 
E^ x (l)=EZx(A) 

= ^( MOiA (AA) + mi.a(AA)) + ^ /i0 . A (AB) = 1 

where 

7 :=E2i(A) =paa (mo,a(AA) + a*i,a(AA)) + ?wo,a(AB) 

denotes the expected total number of parasites in cells of type A in the first generation 
(recall from (|SA1[) that Z = Z (A) = 1). As a consequence, 

n 

E(Z [n] \T [n] = A) = EZ n (A) = /4(1) = J] Eg' Ak (l) = (J)" 

fc=i 

for each neN. It is also well-known that (Z n (A)) n >o dies out a.s., which in terms of 
(Z[„])n>o means that lim„_>. 00 P(Z[ n ] = 0|7]„] = A) = 1, iff 

E log g' A (1) = ^ ( log Mo,a(AA) + log a*i,a(AA)) + ^ log M o.a (AB) < 0. (3) 
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2.2 Properties of #G£(A) and #G 



We proceed to the statement of a number of results on the asymptotic behavior of 
G*(A) and G* conditioned upon Surv(A) and Surv, respectively. It turns out that, if 
the number of parasites tends to infinity, then so does the number of contaminated 
cells. 

Theorem 1. 

(a) J/P(Surv(A)) > and p AA > 0, then P(# G*(A) -> oo|Surv(A)) = 1. 

(b) //P(Surv) > 0, then P(#G^ -4- oo|Surv) = 1. 

Proof. The proof of assertion (jg) is the same as for Theorem 4.1 in [5] and thus omitted. 

(jb| We first note that, given Surv, a contaminated B-cell is eventually created with 
probability one and then spawns a single-type cell process (as EZi(B) > by (|SA5|) V 
Hence the assertion follows again from Theorem 4.1 in if /xb := Mo,B + A«i,B > 1- 

Left with the case /xb < 1, it follows that 



for otherwise, given Surv, only B-parasites would eventually be left with positive prob- 
ability which however would die out almost surely. Next, paa > leads back to jlj) so 
that it remains to consider the situation when j»aa = 0. In this case there is a single line 
of A-cells, namely — > — > 00 — > and (Z„(A))„>o is an ordinary Galton- Watson 
branching process tending P(-| Surv(A))-a.s. to infinity. For n, k £ N, let 



denote the number of B-parasites at generation k sitting in cells of the subpopulation 
stemming from the cell 0™1, where 0™ := 0...0 (n-times). Using pab = 1 and (|SA5[) . 
notably /Ui j a(AB) > 0, /io.B > and /i^b > 0, it is readily seen that 



P(Surv(A)| Surv) = 1, 




Z, 



V 



t)eG„ + fc+ i(B):t)|n+l=O rl l 




and thus 




for all K £ N and k < K. Consequently, 




< PI lim max Z^{n — k, B) = 0| Surv 



\^n->-cxD 0<k<K 




K 



< 




P lim Z K {n - k, B) = 0| Surv(A) 



fc=0 



= 



for all K £ N 



□ 
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The next result provides us with the geometric rate at which the number of con- 
taminated cells tends to infinity. 

Theorem 2. (v^ n =ff= G* (A)) n >o is a non-negative supermartingale and therefore a.s. 
convergent to a random variable L(A) as n — > oo. Furthermore, 

(a) L(A) = a.s. iff E log g' Ai (1) < or v < 1 

(b) P(L(A) = 0) < 1 implies {L(A) = 0} = Ext (A) a.s. 

Proof. That (^~™# G* (A))„>o forms a supermartingale follows by an easy calculation 
and therefore a.s. convergence to an integrable random variable £(A) is ensured. This 
supermartingale is even uniformaly integrable in the case v > 1, which follows because 
the obvious majorant {v~ n # G„(A))„> is a normalized Galton- Watson branching- 
process having a reproduction law with finite variance and is thus L 2 -bounded (see 
Section 1.6 in [4]). Consequently, (v~ n # G* (A)) n >o is uniformaly integrable and 

EL(A) = lim E #G "^ = lim W(Z n {A) > 0), (4) 
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the last equality following from ([2]) in Proposition [2 

As for jig), L(A) = a.s. occurs iff either v < 1, in which case # G* (A) < # G„(A) = 
eventually, or v > 1 and Elog g^ (1) < 0, in which case almost certain extinction of 
(Z„(A)) n >o in combination with (j4]) yields the conclusion. 

(0 Defining r„ = inf{m G N : # G^(A) > n}, we find that 



P(L(A) = 0) < P(X(A) = 0|r n < oo) + P(r„ = oo) 

(A) 



< P fl {#GL fe (A)/^ m ^0} 



fe=l 



r n < oo ) + P(r n = oo) 



< P(L(A) = 0)™ + P(r n = oo) 



for all n > 1, where the #GJ„ ^(A), fe > 1, are independent copies of #G^(A). Since 
P(L(A) = 0) < 1, Theorem □implies 

P(L(A) = 0) < lim P(r„ = oo) = P [ sup#G;(A) < oo ) = P(Ext(A)) 



which in combination with Ext(A) C {L(A) = 0} a.s. proves the assertion. □ 

Since v < 2 and (y~ n j^ G„(A)) n >o is a nonnegative, a.s. convergent martingale, we 
see that 2-"#G*(A) < 2-"#G n (A) -> a.s. and therefore 

^ * #^!), MB -oo. 
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that is, the asymptotic proportion of all contaminated cells is the same as the asymp- 
totic proportion of contaminated B-cells. Note also that 

P(T [W] =A)=E( #C g (A) ) ^0, as,woo. (5) 

Further information is provided by the next result. 

Theorem 3. There exists a r.v. L G [0, 1] such that #G* /2™ — > L a.s. Furthermore, 

(a) L = a.s. iff /Uo,B/Ui,B < 1- 

(ty J/P(L = 0) < 1, then {L = 0} = Ext a.s. 

Proof. The existence of L follows because 2~™#G* is obviously decreasing. As for 
(jllj), suppose first that /xo,bMi,b < 1 and note that this is equivalent to almost sure 
extinction of a random B-cell line, i.e. 

lim F(Z M > O\Z = k,T m = B) = 

n— >oc 

for any k G N. This follows because, starting from a B-cell, we are in the one-type 
model studied in [S]- There it is stated that (Z[ n ])n>o forms a BPRE which dies out 
a.s. iff /Lto,BA*i,B < 1 (see [5j Prop. 2.1]). Fix any e > and choose m G N so large 
that P(T[ m j = A) < £, which is possible by ([5]). Then, by the monotone convergence 
theorem, we find that for sufficiently large K G N 

EL = lim P(Z [n+m] > 0) 

71— tOO 

< lim F{Z [n+m] >0,T [m] =B)+e 

^ J 1 ™, P ( Z [™+™1 > °. Z lm] = k, T [m] = B) + e 

k>0 
K 

^ J™ E P ( Z M > = fc ' T M = B ) + 2e 

k=0 

< 2e 

and thus EL = 0. For the converse, note that 

= EL 

= lim P(Z [n+1] > 0) 

> ItaP^, > 0,T [X] = B)P(Z W > 0|T [0] = B) 

implies = lim JWOO P(Z[„] > 0|T[ ] = B) and thus /io.B/^i.B < 1 as well. 

The proof of (jb| follows along similar lines as Theorem[2j(b| and is therefore omitted. 

□ 
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2.3 Properties of Z n (A) 

We continue with some results on Z n (A), the number of A-parasites at generation n, 
and point out first that (7~ n Z„(A))„>o constitutes a nonnegative, mean one martin- 
gale which is a.s. convergent to a finite random variable W. In particular, EZ n (A) = 7" 
for all n e N . If EZi(A) 2 < 00, 7 > 1 and 

7 := t,E. g ; i (l) 2 =p AA ( M 2 A (AA)+^ A (AA))+p AB Mo,A(AB) < 7, 

then the martingale is further L 2 -bounded as may be assessed by a straightforward but 
tedious computation. The main difference between a standard Galton- Watson process 
and the A-parasite process (Z„(A))„> is the dependence of the offspring numbers of 
parasites living in the same cell, which (by some elementary calculations) leads to an 
additional term in the recursive formula for the variance, viz. 

Yar(Z n+1 (A)) = 7 2 War (Z n (A)) + 7" Var(^(A)) + Cl v n ft(l) 

for all n > and some finite positive constant c±. Here it should be recalled that 
fn(s) = Es Zn< - A K Consequently, by calculating the second derivative of /„ and using 
7 < 7, we obtain 

for some finite positive constant c^. A combination of this inequality with the above 
recursion for the variance of Z n (A) finally provides us with 

00 

Var( 7 -"£„(A)) < 1 + 7~ 2 ^ 7 _fe (Var(Zi (A)) + cic 2 fc) < 00 

fc=0 

for all n > and thus the L 2 -boundness of (7 _ "Z„(A))„> . 

Recalling that (Z n (A)) n > ancl {2 n ) n >o satisfy the extinction-explosion principle, 
the next theorem gives conditions for almost sure extinction, that is, for P(Ext(A)) = 1 
and P(Ext) = 1. 

Theorem 4. (a) If paa — 0, then 

P(Ext(A)) = 1 iff ^o,a(AB) < 1 or v < 1. 

(b) If paa > 0, then the following statements are equivalent: 

(1) P(Ext(A)) = 1 

(2) E# G;(A) < 1 for allneN 

(3) v<\, or 

i/>1, Elog.g Al (l) <0 and q mf i Eg' Al (l) e < 1. 
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(c) P(Ext) = 1 iff P(Ext(A)) = 1 and (i ,b + Mi,B < 1 

Remark. Let us point out the following useful facts before proceeding to the proof of 
the theorem. We first note that, if Elog<7 A (1) < and E<? A (1) log<7 A (1) < 0, then 
the convexity of 6 (->• Kg' A (l) e implies that 

E<7 Ai (l) = inf Eg' Ai {l) e . 

1 O<0<1 1 

If EZi(A) 2 < oo, Geiger et al. [El Theorems 1.1-1.3] showed that 



»(Z n (A) > 0) ~ cnT* fmf Eg' Al (l) e ^ 



as n — > oo (6) 



for some c € (0, oo), where 

re = if Eg Ai (1) log g' Ai (1) < 0, (strongly subcritical case) 

k = 1/2 if Eg A (1) log <? A (1) = 0, (intermediately subcritical case) 

k = 3/2 if Eg' A (1) log <7 A (1) > 0, (weakly subcritical case) 

A combination of @ and ([5]) provides us with the asymptotic relation 

E#G^(A) ~ cn- K v n ( mf ^^(1)^ as n -> oo, (7) 
in particular (with EZi(A) 2 < oo still being in force) 

inf Eg' Ai (l) e <- if supE#G;(A)<oo. (8) 

O<0<1 V „>x 

Proof, (jlj) If paa = and i/ = pab = 1, each generation possesses exactly one A-cell 
and (JU„(A))„>o thus forms a Galton- Watson branching process with offspring mean 
/zo,a(AB) and positive offspring variance (by (ISA3[1 ). Hence a.s. extinction occurs iff 
/Lto,A(AB) < 1 as claimed. If v < 1, type A cells die out a.s. and so do type A parasites. 

"([El T) ^>([E2 ]l " (by contraposition) We fix m e N such that E(#G^(A)) > 1 and 
consider a supercritical Galton- Watson branching process (S n ) n >o with So = 1 and 
offspring distribution 

F(S 1 = k) = F(#G* m (A) = k), fee N . 

Obviously, 

P(5„>fc)<P(#G; m (A)>fc) 

for all k, n E No, hence 

lim P(#G; m (A) > 0) > lim F(S n > 0) > 0, 
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i.e. A-parasites survive with positive probability. 

"jb2])^(|bl|)" If E# G* (A) < 1 for all n E N, then Fatou's lemma implies 

1 > liminf E#G;(A) > E fliminf #G* (A)) 

giving P(Ext(A)) = 1 by an appeal to Theorem[T] 

"|b3l)^(|bTI),(|b2|)" If v < 1 then E#G,;(A) < E#G„(A) = v n < 1 for all neN. So 
let us consider the situation when 

v > 1, Elogg^ (1) < and inf Eg', (if < - 

is valid. By ©, 

E#G:(A) = i/"P(Z„(A) >0) 

for all n e N. We must distinguish three cases: 

Case A. Eg' Ai (1) log g' A (1) < 0. By what has been pointed out in the above remark, 
we then infer 

l=Eg' Ai (l)= inf Eg' Ai (l) 8 <- 
and thus 7 < 1, which in turn entails 

e#g:(a) <EZ„(A) = 7" < 1 

for all n 6 N as required. 

Case B. E 5 ^(l) log 5^(1) > and EZi(A) 2 < 00. Then, by ©, 

P(Z„(A) > 0) ~ en- 3 ? 2 ^m^Eg'^i l) e ^\ as n — >• 00 

holds true for a suitable constant c <G (0, 00) and therefore 

= Km z/ n P(Z„(A) > 0) = liminf E#G;(A) > E (liminf #G* (A)") , 

implying P(Ext(A)) = 1 by TheoremQ] 

Case C. Eg' A (1) log g' A (1) > and EZi(A) 2 = 00. Using contraposition, suppose 

that sup„ eN E#G*(A) > 1. Fix any vector a = (as U) ) u e{o,i},5e{AA,AB,BB} of distribu- 
tions on No satisfying 

ai u J < P (x{^(A, s) = for x > 1 
and u, s as stated, hence 

a (u) > p( X («) (A)S) = ) and < P (xS(A,s)>n) 
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for each n > 0. Possibly after enlarging the underlying probability space, we can then 
construct a cell division process (Z a}V , T v ) ve ^ coupled with and of the same kind as 
(Z v ,T v ) V £t such that 

*SL(A,s) < 45 (M a.s. 



and ¥(xQ tV (A,s)=x) = c$> 



for each u G {0, 1}, s e {AA, AB, BB}, v G T, fc > 1 and x > 1. To have (Z Q ,„, T 1 ,) 1 , eT 
completely defined, put also 

(41,(b),41,(b)) == (4>),4>)) 

for all w £ T and k > 1. Then Z a>1 , < Z„ a.s. and thus 

E< Ai (l) 9 < E. 9Ai (l) 9 , fle [0,1], (9) 

where Z a ^{f\) and <7 Q ,Ai have the obvious meaning. Since the choice of a has no affect 
on the cell splitting process, we have v a = v > 1, while © ensures 

Eio gff ; iAi (i)<Eio gff ; i (i)<o. (io) 

For N G N let a(AT) = (as"" 1 (A0) ue {o,i},se{AA,AB,BB} ^ e tne vec tor specified by 



(*K(a, : 



if 1 < x < N 
if x > N. 



Then EZ Q ( A r) 1 (A) 2 < oo and we can fix TV G N such that sup„ gN E# <G* (jv), ra ( A ) > 1, 
because # ^^^(A) t # G* (A) as AT — > oo. Then, by what has already been proved 
under Case B in combination with (1!))) . (fTU|) and v a (N) > 1, we infer 

inf E.g Al (l) e > inf Eg' A (l) e > -. 
o<e<i 1 o<e<i 1 1 v 

and thus violation of (Ib2l). 



"Jb2])^((b3|" Suppose E#G;(A) < 1 for all n G N and further > 1 which, by 
©, entails lim^oo P(Z„(A) > 0) = and thus E log 5^(1) < 0. We must show 
that Elogg Ai (l) < and inf < 9 <i Eg' Ai (l) e < v' 1 . But given Elog# Ai (l) < 0, the 
second condition follows from ([5]) if EZi(A) 2 < oo, and by a suitable "a-coupling" as 
described under Case C above if EZi(A) 2 = oo. Hence it remains to rule out that 
Elogp A (1) = 0. Assuming the latter, we find with the help of Jensen's inequality 
that 

inf logE 9 ' A (l) e > inf 0Elogo' A (1) = 

or, equivalently, 

inf Eg', {l) e >1>- 

0<9<l aAlV ' - V 
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(which implies info<e<i Kg' A (if = 1). Use once more a suitable "a-coupling" and fix 
a in such a way that 

1 = inf Eg', (if > inf Eg' Ai (if > - 

o<e<i AlV ; o<e<i 3o ' Al1 ' v 

which implies subcriticality of the associated BPRE (Z a ,n(A)) n >o- By another appeal 
to ([5]), we thus arrive at the contradiction 

su P E#G;(A) > supE#Gl„(A) = oo. 

This completes the proof of (0. 

(jg) Since Ext C Ext (A), we see that P(Ext) = 1 holds iff P(Ext(A)) = 1 and the 
population of B-parasites dies out a.s. as well. But the latter form a Galton- Watson 
branching process with offspring mean /Uo,b+A*i,b once all A-parasites have disappeared 
and hence die out as well iff /l*o,b + A*i,b < 1- D 

Theorem 5. Assuming P(Surv(A)) > and thus particularly 7 > 1, the following 
assertions hold true: 

(a) IfEZ 1 (A) 2 < 00 and 7 < 7, t/ien P(Fy > 0) > and EVK = 1. 

(b) IfF(W = 0) < 1, then Ext(A) = = 0} a.s. 

Proof, (gj) As pointed out at the beginning of this subsection, (Z n (A)/j n ) n is a L 2 - 
bounded martingale and thus uniformly integrable. It therefore converges in L 1 to its 
limit W satisfying EW = 1 as well as E>(W > 0) > 0. 

jb| follows in the same manner as Theorem [210 . □ 

3 Relative proportions of contaminated cells 

We now turn to a statement of our main results that are concerned with the long-run 
behavior of relative proportions of contaminated cells containing a given number of 
parasites, viz. 

#{v e g \z v = k} 
Fk{n) •= 

for k G N and n — > 00, and of the corresponding quantities when restricting to con- 
taminated cells of a given type t, viz. 

#{^6G;(t)|z, = fc} 

FkM = — wm — 

for t G {A, B}. Note that 

F k (n) =F k (n,A) + F k (n,B) . 
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Given survival of type-A parasites, i.e. conditioned upon the event Surv(A), our results, 
devoted to regimes where at least one of A- or B-parasites multiply at a high rate, 
describe the limit behavior of Fk(n, A), #G*(A)/#G* and Ffc(n, B), which depends 
on that of Z n {A) and the BPRE Z n (A) in a crucial way. 

For convenience, we define 

P«,t := H-\Z =z,T =t), z e N, t e {A, B}, 

and use E Zjt for expectation under P Zjt . Recalling that P stands for Pi,a, we put 
P* := P(-| Surv(A)) and, furthermore, 

P* it := P, (t (-| Surv(A)) and P" t = F z<t (-\Z n (A) > 0) 

for z £ N and t £ {A, B}. Convergence in probability with respect to P* is shortly 
expressed as P*-lim . 

Theorem |5] deals with the situation when B-parasites multiply at a high rate, viz. 

/^o,b/^i,b > 1, 

In essence, it asserts that among all contaminated cells in generation n those of type 
B prevail as n — > oo. This may be surprising at first glance because multiplication of 
A-parasites may also be high (or even higher), namely if 

^o,A(AA)P A > 1 , A (AA)f- / io, A (AB)f AB > 1, (SupC) 

i.e., if the BPRE (Z n (A)) n >o is supercritical. On the other hand, it should be recalled 
that the subpopulation of A-cells grows at rate v < 2 only, whereas the growth rate 
of B-cells is 2. Hence, prevalence of B-cells in the subpopulation of all contaminated 
cells is observed whenever # G* (B)/# G„(B), the relative proportion of contaminated 
cells within the n th generation of all B-cells, is asymptotically positive as n — > oo. 

Theorem 6. Assuming ^bA'i.b > 1 7 the following assertions hold true: 

(a) 

(b) Conditioned upon survival of A- cells, Fk(n, B) converges to in probability for any 
k G N, i.e. 

P*-lim F k (n, B) =0. 

n— f oo 

Properties attributed to a high multiplication rate of A-parasites are given in Theo- 
rem [71 First of all, contaminated B-cells still prevail in the long-run because, roughly 
speaking, highly infected A-cells eventually produce highly infected B-cells whose off- 
spring m generations later for any fixed m are all contaminated (thus 2™ in num- 
ber). Furthermore, the Fk(n,A) behave as described in [5] for the single-type case: as 
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n — > oo, the number of parasites in any contaminated A-cell in generation n tends to 
infinity and Fk(n, A) to in probability. Finally, if we additionally assume that type-B 
parasites multiply faster than type-A parasites, i.e. 

[1 B ■= Mo.b + Mi,B > 7, 

then type-B parasites become predominant and -Ffc(n, B) behaves again in Bansaye's 
single- type model [5]. 

Theorem 7. Assuming ([SupC]), the following assertions hold true: 

(a) Conditioned upon survival of IK- cells, F k {n, A) converges to in probability for any 
k e N, i.e. 

P*-lim F k (n,A) = 0. 

n— ¥oo 

(b) 

( c) If E^b-Z^ < oo, ^ B > 7 and ^ ,b log ^o,b + Mi,b log A*i,b < 0, then 

P*-lim F k (n, B) = P(^(B) = k) 

n— >oo 

for all keN, where P(^(B) = k) = lim,^^ Pi, B (Z [n] = k\Z [n] > 0). 

4 Proofs 

Proof of Theorem[6} (jg) By Theorem^ 2""#G* -)• L P*-a.s. and P*(i > 0) = 1, 
while Theorem [5] shows that z/~ n #G*(A) — > L(A) P-a.s. for an a.s. finite random 
variable L(A). Consequently, 

#g;(a) _ ( vy ( 2" ^ (#&yj<)\ _ i (V y. #g;(a) ; Q r _ aiS< 



#g: V27 v#g;; V ^ 7 ^ v 2 

as n — )■ oo, for ^ < 2. 

(|b| Fix arbitrary e, (5 > and if S N and define 

( K \ 



D n := <^ F *( n > B ) > 6 ( n Surv(A). 



,fc=i 



By another appeal to Theorem|31 # G* (B) > 2 n L P*-a.s. for all n £ N and L as above. 
It follows that 

#{« € G„(B) : < Z v < K} > 5#G* n (B)l Dn > S2 n Ll Dn , 
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and by taking the expectation, we obtain for m < n 



5E(L1 D J < ±;eI l{o<z e <jf,T„=B}J 

^ ^T E ( S W^B} + #{"eG„:T„| m = A,r„ = B} 

Wg„ / 
^ 4 E P (° < K ' T *\™ = B ) + i E#G m (A) 



on V " — --i-Ti"' -/ ' or. 

•uSG„ 



■DC 



< y^ / y> P(Z M = z,T M = B) \ / P,,b(0<^ <-g) \ /yy 

- I 2^ 2 m / I ^ 2™-™ / V2/ 



z =\ \uei 

oo 



< ^2nZ [m] = z, T [m] = B) P,, B (0 < Z [n _ m] < K) 



Since v < 2 we can fix m € N such that (i//2) m < e. Also fixzo€N such that 

P(Z [m] > z ) < e. 

Then 

5 E (L 1 D „ ) < J2 P ( Z M = T M = B) P,,b (0 < %-m] < ^) - 



2>1 



2=1 

But the last sum converges to zero as n — ► oo because, under P Zj b, (■^[«])«>o is a 
single- type BPRE (see [5]) and thus satisfies the extinction-explosion principle. So 
we have shown that EL ljj n — > implying P(D n ) — > because L > on Surv. This 
completes the proof of the theorem. □ 

Turning to the proof of Theorem [7J we first note that part (a) can be directly inferred 
from Theorem 5.1 in [3] after some minor modifications owing to the fact that A-cells 
do not form a binary tree here but rather a Galton- Watson subtree of it. Thus left 
with the proof of parts (b) and (c) , we first give an auxiliary lemma after the following 
notation: For v € G n and k £ N, let 

G* k (t,v) := {u G G* n+k (t) : v < u} 

denote the set of all infected t-cells in generation n + k stemming from v. Let further 
be 

G„(A,B) :={«6C; +1 (B) :T u| „=A} 
which is the set of all infected B-cells in generation n + 1 with mother cells of type A. 
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Lemma 1. If ( |SupCD holds true, then 

P,„» *M = „ > 0, 

™^°° #g*(a) 

w/iere /3 := lim^oo E ZiA #G*(B). 

Proof. Since z h-> E Z! a#G^(B) is increasing and Ei j a#Gj(B) > by our standing 
assumption (|SA5I) . we see that (3 must be positive. Next observe that, for each neN, 

#g;(a,b)= # g i( b ^)> 

ueG*_ 1 (A) 

where the # GJ(B, i>) are conditionally independent given Z n (A) > 0. Since # G* (A) — > 
oo P*-a.s. (TheoremHJ) and P" = P(-|Z„(A) > 0) A P*, it is not difficult to infer with 
the help of the SLLN that 

7r „v ; uSG*(A) 

P* 

where a n ~ b n means that P*(a n /b n — > 1) = 1. Now use E Zj a#GJ(B) f /3 to infer the 
existence of a zq S N such that for all z > zo 

E z ,a#G*(B) > 0(1 — e). 

After these observations we finally obtain by an appeal to Theorem (7{a) that 

e * E E ^^#^(B) 

2 >^o rrl nK J\ v - u; „ e { tte6 * (A)|Z„>«q} 

> 0(1 -e) £ ^V(n,A) 

2>z 

— > 0(1 - e), n — > oo. 
This completes the proof of the lemma. □ 

Proof of Theorem [7](b) and (c): Let e > and JVeN. Then 

n-l 

#g;(b) = E E #G*_ fc -i(B,«) 

fc=0 t,eG*(A,B) 
n-l 

> E E #G;_ fe _ 1 (B, v ) 

fe=0 «G{«GG*(A,B)|Z„>z} 
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#{v € G;(A, B):Z v > z } 

#g;(a,b) 

This together with Lemma Q] yields 

#{v€G* n (A,B):Z v >z } 

#g;(a) 

and thereupon 



1. 



(11) 



> ^ #G^(B, V ) 

«e{ueG*_ 1 _ m (A,B)|z u >z} 

a.s. for all n > to > 1 and 2 £ N, and thus 

> = (Jv#g;(a) > #g;(b)) 

< p* U#g;(a)> £ #c; n (B,«) 

\ w6{ueG*_ 1 _ ro (A,B)|Z„>z} 

Fix to so large that 

4iV 

2 m (l -£)>—• 

Then, since 

lim P Z)B (#G^ = 2 m ) = 1, 

z— >oo 

there exists zq £ N such that 

F,,b(#G; u = 2 m ) > 1-e 

and therefore 

E 2 , B #G; n > (l-e)2 m > (12) 

p* 

for all z > zq. Moreover, ^2 k>z Ffc(n, A) 1 by part (a), whence 



„^eG;(A,B):Z„> Zo y, f (i:s) 



#g;(a) - 2, 

for all n > no and some no G N. By combining (jlip and (|13[) . we now infer for all 
71 > no + to 

#g;(a) > i 



#G* AT + 1 



< P* (tV#G:(A)> £ #G^(B,«) 

uS{«GG;_ 1 _ m (A,B):Z u >z} 
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2N S w e{ ue G;_ 1 _ TO (A,B):Z u > 2; } # G m( B > w )\ 
T > #{u G G;_ 1 _ m (A, B) : Z M > z} J + 6 

„ / 2iV E 1 , G { Me G;_ 1 _ m (A,B):Z u > z } # G m( B , v ) \ P(Z„_ m (A) > 0) 
\T > #{u i G;_ 1 _ m (A, B):Z U > z} ) P(Surv(A)) 

„ 2N E t { r G: — (A|B)!Z ^' } ft, m (zo) \ F(Z„_ m (A)>0) 
^ /3 #{« G G^.JA, B) : Z u > z} ) P(Surv(A)) 



where the Si,m(^o) are iid. with the same law as #{i> G G* rl (B) : Z = zo,?0 = B}. 
The LLN provides us with n\ > uq + m such that for all n > ri\ 



/^#{«eG;_ 1 _ m (A,B):Z u > 2 } 



>Eg i , m (z )/2 > 1-e. 



y #{it G G*_ x _ m (A, B) : Z u > z} 
By combining this with (|12j) . we can further estimate in the above inequality 



#<g;(a) > i 



#G* JV + 1 

I — > E&,m(«>)/2 > — j + £ j P(Sur v(A)) + 6 



/ P(Z n _ m (A) > 0) \ „_ 
^ P(Surv(A)) ) e 



J -> 2e. 



This completes the proof of part (b). 

As for (c), we will show that all conditions needed by Bansaye [3] to prove his 
Theorem 5.2 are fulfilled. Our assertions then follow along the same arguments as 
provided there. 

Step 1: (/ig"Z„(B)) n >o is a submartingale and converges a.s. to a finite limit W(B). 
The submartingale property follows from 



E(Z„ +1 (B)|Z„(B)) = E 1^2 ( 

Z V Q 1{T„ = B} +Z v i 1{T„ 1 = B} ) 



i vet. 



Z n (B) 



= Z n (B)E[x(°)(B)+X( 1 )(B)j +E (^o1{t„„=b}+2,i1{t„ 1 =b } ) 

\ueG*(A) 

> Zn(B)nB 

for any neN, while the a.s. convergence is a consequence of 

,'ZJB) 
sup E | n \ ' | < oo 



Z n (B) 



ign \ Mb 
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which, using our assumption 7 < /zb, follows from 



-n+H 



Mb / \ Mb 



\ Mb «eg»(A) / 



Mb / Mb +1 



Z„(A)E(Z l-{To = B} +^1 1{T 1 = B}) 



Z„(B)\ ^ab / 7 V 



E 

Mb J Mb \Mb, 

n / \ k 00 / \ k 

Mab / J_\ < Mab / 7 \ 

^ ^ W mb toW 



< 00 



for any n£N. 



Step 2: {W(B) = 0} = Ext a.s. 
The inclusion D being trivial, we must only show that P(W(B) > 0) > P(Surv). For 
i > 1, let (Z i: „(B)) n > be iid copies of (Z„(B))„> under P ljB . Each (Z i>n (B)) n > 
forms a Galton- Watson process which dies out iff /Xg n Zj )n (B) — > (see [5]). Then for 
all m,N GN, we obtain 

F(W(B) > 0) = P ( lim ^ m +" ( „ B) > 
V n ^°° Mb 

v J I Inn — V >() 



lim N 

\"^°°Mb ^ 

/ _ 1 ^(B! 

lim } 

[ Um y > 0, Z m (B) > iV ] 

Mb / 



Mb 

> P|Jim -"y: B) ^p>0,Z m (B)>7V 



> 



> P (Z m (B) > JV) - Pi, B f lim V = 

= P (Z m (B) > N) — Pi, B f lim ^4r^ = ° N 
= P (Z m (B) > JV| Surv) P(Surv) - P ljB (Ext)* 



and then, upon letting m and iV tend to infinity, 

P(W(B) > 0) > P(Surv) 
because Pi,b (Ext) < 1 and by Theorem [1] 



22 



Step 3: sup n > E£„ < oo, where £ n := (mb/2) n Zy n y 
First, we note that (Zr n i) n >o, when starting with a B-cell hosting one parasite (under 
Pi,b), is a BPRE with mean /zb/2 (see [5]). Second, we have 

EZ [n] l {T[n]=A} = P(T [n] =A)EZ„(A) = (J)™ 

and thus 

n-l 

EZ[„] = EZ[„] 1{t w =a} + ^ EZ[„] l{T [m] =A,r [m+1] =B} 

m— 

n-l 

= (2) + l{T [m] =A} E 1.A^[1] 1{T [1] = B} Ei ; BZ[ n _ m _i] 

m=0 
n-l 

= (£+•>£ 



\ m / ^ig \ n— m— 1 

2 



m=0 

for all neN, where 77 := E^a-Z^j 1{ T B |. This implies 

n-l 



+?E £ <» (») 



supE£„ = 

neN \VB/ MB^VMb/ ^ VMB 

for some c < 00. 

Step 4: lim^oo sup„> E£„ l{ Z[n] >K} =0. 
By our assumptions, (Zr n i)„>o, when starting in a B-cell with one parasite, is a strongly 
subcritical BPRE with mean /x B /2 (see [S]). Hence, by [I] Corollary 2.3], 

Jim supEi jB £n 1{z m >k} = 0, (15) 

K^oc n>0 



which together with ([14")) implies for n, m 6 N 
Jim supE£ n+TO l {z , n+m]> K] 

Jf-KX>„> 

< lim supE^„ +m l {Z[n+m]>Rr }l{ T[m]=B }+supE^ n+m l {T[m]=A } 

K— >oo n >Q n>0 

< lim SUpEi B ^n 1{Z,„,>K} E^ m +E^ m 1/ T =A } supE^„ 

< (-L) SU P E£„ 
VA*B/ neN 

and the last expression can be made arbitrarily small by choosing m sufficiently large, 
for 7 < This proves Step 4. 

Final step: Having verified all conditions needed for the proof of Theorem 5.2 
in [5] , one can essentially follow his arguments to prove Theorem [3(c) . We refrain 
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from supplying all details here and restrict ourselves to an outline of the main ideas. 
First use what has been shown as Step 1 - 4 to prove an analogue of [5j Lemma 6.5], 
i.e. (control of filled-in cells) 

f #{v 6 G* + (B) : Z v \ n > K} \ 
hm sup P q > rj = (16) 

^->°°n,g>0 V #G n+g (B) / 

for all r\ > 0, and of Prop. 6.4], i.e. (separation of descendants of parasites) 

v p * ( #{v £ G* n+q (B) : Z„, n < K,N n (v) > 2} \ 

hm supP — — — — >n 1=0 (17) 

^°°«>0 V #G n+(? (B) / 

for all 77 > 0, -K" S N, where N n (v) denotes the number of parasites in cell v\n with 
at least one descendant in cell v. In particular, (|16p (with q = 0) combined with 
# G* (B) -> 00 P*-a.s. implies the existence of a K" > such that for all N e N 



lim inf P* V Z„ 1 {Z . <X} > JV =1. (18) 

I-UGGJ(B) 



Using (fl"6"|) and (jTTJ) , we infer that, for all ry,e > 0, there exist K\ > Kq and qo € N 
such that for all n € N 



F k (n + q ,B) 



#{v € G* n+qo (B)\Zv = k,Z v]n < K u N n (v) = 1} 
#{« G K+ qo ( B )\ Z v\n < K 1: N n (v) = 1} 

V v ' 



> f}\ < e. 



Since # G*(A)/# G£(B) 0, we further get 

p* #{^ £ G; +g0 (B)K = fc,Z„| n < K u T v \ n = B,N n (v) = 1} 



#{v G G* (B)|Z„, n <K u T vln = B,N n (v) = 1} 



asn-> 00, which puts us in the same situation as in the proof of [5J Thm. 5.2]. Now, 
by using (fl~8|) and the LLN, we can identify the limit of J„ which is in fact the same as 
in Step 1 of the proof of [5J Thm. 5.2]. A reproduction of the subsequent arguments 
from there finally establishes the result. □ 
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5 Glossary 



T cell tree 

G„ set of cells in generation n 

G n (t) set of cells of type t in gegeration n 

G* set of contaminated cells in generation n 

G* (t) set of contaminated cells of type t in generation n 

T v type of cell v 

p s probability that the daughter cell of an A-cell is of type s 

po probability that the 1 st daughter cell of an A-cell is of type A 

pi probability that the 2 nd daughter cell of an A-cell is of type A 

v mean number of type-A daughter cells of an A-cell 

(X(°\A, s), X^(A, s)) offspring numbers of an A-parasite with daughter cells of type 

s € {AA, AB, BB} 

(X^(B),X^(B)) offspring numbers of a B-parasite 

Z v number of parasites in cell v 

jUj,t(s) mean number of offspring of a t-parasite which goes in daughter 

cell i e {0, 1} if daughter cells are of type s <E {AA, AB, BB} 

yUi.B mean offspring number of B-parasites which go in daughter cell 

ie{o,i} 

reproduction mean of a parasite in a B-cell 
Z n number of parasites in generation n 

Z n (t) number of parasites in t-cells in generation n 

Ext / Surv event of extinction/survival of parasites 

Ext(t)/ Surv(t) event of extinction/survival of t-parasites 

Z[ n ] number of parasites in a random cell in generation n 

Z n (A) number of parasites of a random A-cell in generation n 

f n (s\A), f n (s) quenched and annealed generating function of Z n (A), respec- 

tively 

g\ n (s) generating function giving the n-th reproduction law of the pro- 

cess of a random A cell line 
7 mean number of offspring of an A-parasite which go in an A-cell 

7 := isEg' Ai (l) 2 = p AA >A (AA) + /i?, A (AA)) + Pab Mo,a( AB ) 

P Zjt probability measure under which the process starts with one t- 

cell containing z parasites 

the same as before but conditioned upon Surv (A) 
the same as before but conditioned upon survival of A-parasites 
in generation n 



■ z,t 
z.t 
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